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In this work we study the QQnn systems from both four-quark and meson-meson molecule for- 
malisms, arriving to consistent conclusions when the proper bases are chosen. These states are 
interesting as they are thought to be found in experiments soon at LHC or RHIC. Such a source 
of experimental information will significantly improve our current understanding of charm spec- 
troscopy. 
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After the discovery of in 1974 [jl|, three quiet decades followed in charmonium spec- 
troscopy. Although many states belonging to the cc spectrum were discovered, all of them nicely 
fitted into the slots theoretically predicted. The situation was analogous for D and D s mesons: their 
spectra filled up without difficulties. However, two new discoveries performed in 2003 meant a 
turning point. The first one, a charm-strange state D* 7 (2317), was reported by BABAR ^ and 
promptly confirmed by different collaborations. Its small width was not compatible for such a 
small mass if assumed a standard qq configuration. The second one is the celebrated X(3872), first 
reported by BELLE [|J]. Although both states are firmly stablished, they can hardly be understood 
in a simple qq scheme, what makes their nature not being definitively settled. 

By having a look at the charmonium spectrum, one can check that things look easy below the 
threshold for the DD production: all the states match the predictions for the cc spectrum based on 
a simple color Fermi-Breit interaction. Above this threshold, it gets all more involved, with many 
puzzles awaiting solution. There are still empty levels of the cc spectrum for whom a candidate 
has not been found, whereas several of the abundant new states do not fit in the cc spectrum, as the 
confirmed X(3872) or Y(4260), or the still dubious Z(4050) or Z(4250). These latter are of par- 
ticular interest since they have been interpreted as charged resonances and therefore cannot belong 
to the charmonium spectrum. Thus we have a zoo of particles named X, Y or Z whose structure 
has yet to be determined. Some possible explanations have been postulated so far: molecular, cc + 
molecular, hybrid, or cc + higher order Fock space components, among others. The wave function 
of a meson (B = 0) within the constituent quark model, where explicit gluon degrees of freedom 
are frozen in terms of a quark constituent mass, can be written as: 

|¥(fl = 0)) = OCi\qq) + a 2 \qqqq) + ... (1) 

where |a;| 2 = 1. Ground state qq mesons have negative parity. Positive parity mesons can be 
reached either through a four-quark configuration in a L = state or by adding a unit of orbital 
angular momentum to the qq pair. In the constituent quark model, the mass of such a pair is around 
600 MeV []|], whereas the mass shift produced by a unit of L stands around 500-600 MeV [g]. 
None of these two mechanisms is suppressed by the other and therefore the four quark piece is not 
negligible now. The challenging consequence of this so-called unquenching of the quark model 
[J5|], that seems to be unavoidable nowadays when pursuing a description of the excited hadronic 
spectrum, is the appearance of exotic states. They could be stable in nature, with a mass that is 
similar to the positive parity excitations, and they could not be described by the lowest order Fock 
space components of the naive quark model, \qq). This would be the case of doubly charmed or 
bottom mesons ( QQnn states, where Q stands for the heavy quark and n for a light one), long ago 
suggested in the literature [[§]. They were found to be stable against dissociation provided the ratio 
between the heavy and light masses was large enough (in a clear contrast with QQnn, for instance). 
If such states were below the DD threshold, they would be narrow and would clearly show up in 
the experiments. 

Let us get more insight into these QQnn systems by approaching them from two different 
techniques. We solved in first place the Schrodinger equation through the Hyperspherical Harmonic 
(HH) method [St]. This is done by expanding the trial wave function in terms of HH functions, 
generalizing thus the well-known Spherical Harmonic formalism. The main difficulty one has to 
fight with is the construction of base states with the proper symmetry. The Pauli principle needs 
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to be imposed since some of the quarks in QQnn are identical. A constituent quark cluster model 
(CQC) has been employed for the interaction. When screening all the J p channels for the cchh 
system, only one state, the (I)J P = (0)1 + , is found to lie below the DD threshold. 

To get a deeper understanding on the structure of this four-quark state, one can develop a 
method to compute probabilities of physical vectors in any four-body system [[l0|]. It has been 
found that exotic states would present significant probabilities of different physical channels. There 



are three different ways to couple qqqq to a colorless state: 

[(qiqi) (9394)] = {!3i23 3 4),|6 12 6 3 4)} = { |33)' 2 , |66)' 2 } , (2) 

l(qm) to)] = {[I13I24), [8 13 8 24 )} = {|ll) c ,|88) c }, (3) 

\{ qi q A ) (q 2 q 3 )} = {|1 14 1 23 ), |8 14 8 23 )} = {|l'l') c , |8'8') c }. (4) 



Bases (Q) and (Q) contain physical vectors, i.e. the [ 1 1) and | l'F) components, together with octet- 
octet (hidden color) ones. All vectors inside each basis are orthogonal, however bases are not 
orthogonal to each other. We aim to express the (0) 1 + four-quark state in terms of the physical 
components. When solving the problem in base ([?]) we found a probability of 62.6% to obtain a 
hidden-color nonphysical vector |88) c . However this piece can be decomposed again, as bases (||) 
and (Q) can be connected through the following antiunitary transformation: 

|ll) c = cosa|l'l') c + sina|8'8') c (5) 
|88) c = sina|l'l') c -cosa|8'8') c . (6) 

That would generate another physical state, | l'l'^, together with anew hidden-color piece (|8'8') c ). 
This latter can be reexpanded again in terms of the first base through Eq. @ and this process 
should be repeated up to infinity getting a series that can be summed up JTo|] - After doing so one 
gets 50.5% of DD* and 49.5% of D*D*, in our present (I)J P = (0)1+ case, being |ll) c = \DD*) 
and |l'l') c = \D*D*). One may compare the probabilities obtained for the only bound state of the 
QQnn system with those for an unbound channel, the (I)J P = (1)1 + for instance. We find there 
is no probability associated to the D*D* channel. This comparison helps us to see the relation be- 
tween the binding energy and the compactness of the system. The higher the binding energy, the 
more compact gets a system and the more spread the probabilities over different physical channels 
are. 

Let us now face the problem from the two-meson interaction point of view. For that purpose 
we start from a physical system made of two D mesons, M\ and M 2 (Mi = D,D*), with quantum 
numbers (I)J P in a relative S state. They interact through a potential V that contains a tensor 
force. Then, in general, there is a coupling to the M\M 2 D— wave and to any other two D-meson 
system (DD, DD*, D*D*) that can couple to the same quantum numbers (I)J P . Thus, if we denote 
Di = DD, D 2 = DD* and D 3 = D*D*, the Lippmann-Schwinger equation for the M\M 2 scattering 
becomes 

tap-Ji (P<*iPP> E ) = V aP;ji (P«>PP)+ L L L Py d Py V a T ,ji (P«,Py) 

y=D k ^=°> 2 ° 

(£=1,2,3) 

x G r (E;p r )/ 7 r ^ pSp (p r ,pp-,E) , a,P=D h D 2 ,D 3 , (7) 
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Figure 1: Fredholm determinant for the (I)J P = (0)1 + ccnn system. The dashed line corresponds to the 
calculation considering only DD* singlet-singlet color states. The dashed-doted line stands for the results 
considering only the D*D* system. The solid line represents the results including the full set of color states 
or, in other words, the DD* — D*D* coupled system (see text for details). 



where t is the two-body scattering amplitude, j, i, and E are the angular momentum, isospin and 
energy of the system, l a Sa> ^ySy, and Ipsp are the initial, intermediate, and final orbital angular 
momentum and spin, respectively, and p r is the relative momentum of the two-body system y. The 
basic ingredient to solve the scattering problem are the interacting potentials. They are taken from 
the CQC model, the very same that we used for the Hyperspherical Harmonic formalism. To obtain 
them from the basic qq interaction we use a Born-Oppenheimer approximation. 

All positive-parity channels made by S-wave interacting mesons, up to J p = 2 + , have been 
analyzed. Again, only one of them is attractive enough to be bound. This channel is (I)J P = (0)1 + , 
the same that was found to bind from our four-quark study. We can have a look at the potentials 
[11] and check that both DD* and D*D* are attractive. It gets all more quantitative and clear 
when attending at the Fredholm determinant, as done in Fig. [T]. Although both total potentials, 
DD* and D*D* are attractive (Fredholm determinant smaller than 1), neither of them is bound. 
To achieve a bound system, or a Fredholm determinant equal to zero, one needs to consider the 
transition between DD* and D*D*, and solve the coupled-channel system. Again the coupling, 
as we saw from the HH formalism is the main ingredient to make the system bound. At short 
distances, identical quarks can recouple to different vectors of the Hilbert space. If any of them 
is not considered, a fundamental ingredient of the calculation would be neglected. Such an effect 
would never happen when dealing with hadronic degrees of freedom. In that case, we would just 
have an additional channel 275 MeV above, not giving any significant contribution in a coupled- 
channel calculation. 

The need to incorporate the complete Hilbert space is also evident when the four-quark for- 
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malism is employed: if the higher D*D* channel was not included in the calculation, one would be 
artificially neglecting the large portion of vector |88) c in base (||). If the model parameters were fit- 
ted to some observables they would necessarily include the effect of such a restricted Hilbert space 
and one would never be able to use them for different quantum numbers without arriving to wrong 
conclusions, loosing thus any predictive power. We can therefore conclude that formalisms based 
on four-quark and meson-meson configurations are fully compatible whenever all the relevant 
basis vectors are taken into account. 

When moving to the bottom sector one finds four more canditates for observation. They are 
(I)j p = (0)1+, (0) 0+, (1) 3" and (0) 1~. Together with the (0)1+, the interest in these states is 
increasing as for the first time there are chances to observe such large mass exotic states in quite 
a near future: LHC may discover tetraquark states via gluon-gluon fusion due to both large num- 



ber of events and their unique signature in the detectors [|12j]. Also at RHIC the identification of 



hadronic molecular states by means of relativistic heavy ion collisions has been suggested by em- 



ploying the coalescence model for hadron production [ ]13| , |14| ]. If any of these suggestions gets real, 
the new data related to these double charm or bottom exotic systems will give a huge contribution 
to the heavy quark spectroscopy. 
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EU FEDER under Contract No. FPA20 10-2 1750, and by the Spanish Consolider-Ingenio 2010 
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